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^ ABSTRACT 

o 

j—j We report sub-arcsecond scale radio continuum observations of a field of 35 by 

22 arcmin centred in NGC 3351 obtained with the Multi-Element Radio Linked Inter- 
ferometer Network (MERLIN). We found 23 radio sources in this field, 6 of which are 
projected within the D25 isophote of the galaxy, and 3 are located inside the central 

y-H 100 arcsec in radius. Two of these three are significantly extended, while the third one 

^ is relatively compact. This one is the only source with a previously detected coun- 

terpart at other wavelengths and could constitute the radio counterpart of a young 

ON supernova remnant. The other two are probably related to jets from a background 

AGN. We are not able to detect individual supernovae or SNRs in the central region 
(r < 600 pc) of the galaxy. This could imply that the ionising populations of the 

^j- circumnuclear star-forming regions (CNSFRs) are too young (less than a few Myr) 

£^ to host supernovae. Also the diffusion length of the relativistic electrons in the ISM 

associated with the SN from previous events of star formation seems to be larger than 
our maximum resolution of 50 pc in radius. Detecting the thermal bremsstrahlung 
emission from the circumnuclear Hll regions probably requires deeper observations. 

k> Key words: Galaxies: individual: NGC 3351 - radio continuum: galaxies - galaxies: 

ij starburst. 



1 INTRODUCTION 

The radio emission of a normal spiral galaxy consists of a 
thermal bremsstrahlung component, arising from inelastic 
Coulomb collisions between the free electrons and ions of the 
interstellar medium (ISM), and a synchrotron component 
emitted by relativistic cosmic-ray electrons from supernova 
remnants (SNRs) and active galactic nuclei (AGN) spiralling 
in the galactic magnetic field (see e.g. Condon 1992). Star- 
forming regions may host a large number of supernovae, and 
the Hn regions created by young massive stars can be de- 
tected as thermal sources (e.g. Kobulnicky & Johnson 1999; 
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Turner et al. 2000; Johnson et al. 2001, 2003, 2009; Can- 
non & Skillman 2004; Tsai et al. 2006, 2009; Rosa-Gonzalez 
et al. 2007). In particular, the inner (~ 1 kpc) parts of some 
spiral galaxies show higher star formation rates (SFRs) than 
usual, frequently arranged in a ring pattern around the nu- 
cleus. Circumnuclear star-forming regions (CNSFRs), also 
referred to as "hot-spots" (Sersic & Pastoriza 1965, 1967), 
are similar to luminous and large disk Hn regions, but they 
are more compact and show higher peak surface brightness 
(Kennicutt et al. 1989). Their massive stars can dominate 
the observed ultraviolet (UV) emission even in the presence 
of AGN (Gonzalez Delgado et al. 1998; Colina et al. 2002), 
and the ionised gas can contribute up to 50 per cent to the 
total H/3 luminosity in some Seyfert 2 nuclei (Cid Fernan- 
des et al. 2001). The Ha luminosities, typically higher than 
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10 39 erg s 1 , overlap with those of Hn galaxies (see e.g. Mel- 
nick et al. 1988; Diaz et al. 2000; Hoyos & Diaz 2006, and 
references therein). Due to their location near the centre of 
the galaxy, CNSFRs display high metal abundance (see e.g. 
Diaz et al. 2007, and references therein), which makes them 
ideal laboratories to investigate star formation in environ- 
ments of large metallicity. 

NGC3351, also known as M95 and UGC5850, is a 
nearby (10.05 Mpc, 49 pc/arcsec, Graham et al. 1997) 
SBb(r)II early-type barred spiral galaxy (Sandage & Tam- 
mann 1987) in the Leo group. Its coordinates are CK2000 = 
10 h 43 m 57!7, feooo = +11° 42' 12 / /7, and its major- and 
minor-axis diameters at the 25 magnitude isophote (D25) 
are 7.4 x 5.0 arcmin, with a position angle (PA) of 13° 
(RC3 catalogue, de Vaucouleurs et al. 1991). NGC3351 has 
a star-forming circumnuclear ring, and has been classified 
as a "hot-spot" galaxy by Sersic & Pastoriza (1967). Early 
detailed studies of its nuclear regions by Alloin & Nieto 
(1982) showed that NGC3351 harbours high-mass star for- 
mation along a ring of about 20 arcsec in diameter in its cen- 
tral zone. In fact, it can be considered a nuclear starburst 
galaxy, since the SFR per unit area is significantly higher in 
the nuclear region than in the disc (Devereux et al. 1992). 
HST UV images show that the present star formation in 
the ring of NGC 3351 is arranged in complexes of diame- 
ters between 1.4 and 2.0 arcsec, made up of several high- 
surface brightness knots a few parsec in size embedded in a 
more diffuse component (Colina et al. 1997). A circumnu- 
clear SFR of 0.38 M yr _1 was inferred by Elmegreen et al. 
(1997) from near infrared (IR) photometry in the J and K 
bands, and a value of 0.24 yr _1 was derived by Plane- 
sas et al. (1997) from the Ha emission. From CO observa- 
tions, these authors estimated a mass of molecular gas of 
3.5 x 10 8 M© inside a circle of 1.4 kpc in diameter. Accord- 
ing to Jogee et al. (2005), NGC 3351 hosts 5.3 x 1O 8 M of 
molecular hydrogen, with most of the emission coming from 
the inner ~ 600 pc. Within this radius, a SFR of 0.5 M© yr _1 
can be estimated from the Br7 luminosity given by Puxley 
et al. (1990), whereas the non-thermal component of the ra- 
dio continuum emission favours a lower value of 0.3 M© yr _1 
(Jogee et al. 2005). van de Ven & Chang (2009) estimated 
that the molecular mass inside the circumnuclear ring is 
about 5 x 1O 7 M , and Leroy et al. (2009) found a well- 
defined molecular ring with a diameter of about 50 arcsec, 
located at the end of the stellar bar (with a semimajor axis 
of 47 arcsec; see e.g. Martin 1995). These CO structures are 
compatible with the excesses observed by Regan et al. (2001) 
and Heifer et al. (2003). A small CO bar perpendicular to 
the large-scale stellar bar was also observed by Devereux 
et al. (1992) and Heifer et al. (2003). 

NGC 3351 is part of the SINGS Legacy Science Program 
sample (Kennicutt et al. 2003), and it has been included in a 
wide variety of statistical studies. The SINGS data set pro- 
vides broadband imaging, as well as low- and high-resolution 
spectral maps in the IR. The integrated spectral energy dis- 
tribution (SED) of NGC 3351 from 0.15 to 850 /xm has been 
derived by Dale et al. (2007) by combining the Spitzer data 
with those from 2MASS and IRAS in the IR, GALEX in 
the UV, and other data in the optical, modelling the dust 
and stellar contributions. They concluded that most of the 
sample's spectral variations stem from two underlying com- 
ponents, one representative of a galaxy with a low IR-to-UV 



ratio and one representative of a galaxy with a high IR-to- 
UV ratio, as given by a principal component analysis of the 
sample. Draine et al. (2007) used the IRAC and MIPS ob- 
servations to estimate the total dust mass (1O 7-46 M0) and 
the fraction of the dust mass contributed by PAHs (3.2 per 
cent). They found that the dust properties of spiral galaxies 
resemble those in the local environment of the Milky Way, 
with similar dust-to-gas ratio and similar PAH abundance. 
Munoz-Mateos et al. (2009) found that NGC 3351 shows sig- 
nificantly larger IR concentration indices than the remaining 
galaxies of the same morphological type of the SINGS sam- 
ple. 

We have studied the circumnuclear region of NGC 3351 
in two previous works. In Diaz et al. (2007), the physical con- 
ditions of the gas were determined for 7 CNSFRs. We devel- 
oped a semi-empirical method for the derivation of chemical 
abundances in the high-met allicity regime, obtaining values 
consistent with solar within the errors and comparable to 
those found in high-metallicity disc Hn regions using di- 
rect measurements of electron temperatures. An analysis of 
the kinematics of the gas and stars in 5 CNSFRs and the 
nucleus of NGC 3351 is presented in Hagele et al. (2007). 
We derived the dynamical masses (between 4.9 x 10 6 and 
4.3 x 10 7 M©) for the stellar clusters in these regions, and 
found that the CNSFR complexes, with sizes of about 100 
to 150 pc in diameter, are seen to be composed of several 
individual star clusters (with sizes between 1.7 and 4.9 pc 
on a HST image). The radial velocity curve shows deviation 
from circular motions for the ionised hydrogen, consistent 
with infall towards the central regions of the galaxy at a 
velocity of about 25kms _1 (in agreement with Rubin et al. 
1975; Dicaire et al. 2008). These last authors also point out 
that the Ha velocity field is fairly regular outside the central 
bar, and Jogee et al. (2005) found that the velocity field of 
the molecular material in the inner 500 pc is generally dom- 
inated by circular motions, with some evidence of weaker 
non-circular streaming. 

In this paper, we will look for radio-continuum signa- 
tures of massive star formation in NGC 3351, both in the 
nucleus (in particular, the CNSFRs) and the galaxy disk. 
Section 2 describes the observations and the data reduction. 
The spatial distribution of radio sources and their observed 
properties are discussed in Section 3, and Section 4 is de- 
voted to the interpretation of these results in terms of star 
formation activity. Our main conclusions are briefly sum- 
marised in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed NGC 3351 at 1.42 GHz on 2008 Jan 26-28 us- 
ing 6 antennas of the Multi-Element Radio Linked Interfer- 
ometer Network (MERLIN). The Cambridge telescope has 
a diameter of 32 m and provides the longest baselines of 
up to 217 km; the other five antennas have effective diame- 
ters of 25 meters. The field centre was at right ascension 
10 h 43 m 57!7000, declination +ll°42 / 14 / /000 (J2000). We 
used the point-like QSO B 1023+ 131, at a separation of less 
than 5°, as the phase reference source, with a duty cycle of 
about 6.5:1.5 min on the target and phase reference, respec- 
tively. The position used for B1023+131 (10 h 25 m 56^2854, 
+ 12°53 / 49 / /024) was taken from Browne et al. (1998). We 
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Table 1. Sources around NGC3351 detected by MERLIN. The first column gives the adopted reference name for each source (in 
brackets the adopted names for the three closest sources to the centre of NGC3351), the second column gives the FIRST counterpart 
name (within 1 arcsec of the MERLIN source in all cases) and the third column gives the radial offset of the MERLIN sources from the 
pointing centre. The remaining columns give the MERLIN source properties: positions, peak flux densities, integrated flux densities, and 
their corresponding errors. These measurements for Rl, R2 and R3 were obtained by fitting a single component; see Table 2 for further 
details of these sources. 



Source 


FIRST 


Radius 




R.A. 




Dec. 


CRA 


CDec 


Peak 


Int. 


0"rms 




name 


(arcmin) 






(J2000) 






(mas) 


(mas) 


(mJy b- 1 ) 


(mJy) 


(mJy b- 1 ) 


1 


J104252. 6+114853 


17.283 


10 


42 


52.5564 


+11 


48 


54.157 


47 


76 


0.31 


1.78 


0.08 


2 


J104253. 2+114848 


17.076 


10 


42 


53.3084 


+11 


48 


48.482 


27 


44 


0.55 


1.91 


0.08 


3 


J104254. 6+114841 


16.727 


10 


42 


54.6482 


+11 


48 


41.204 


24 


38 


0.63 


5.33 


0.08 


4 


J104301. 1+115230 


17.228 


10 


43 


01.1845 


+11 


52 


30.364 


27 


44 


0.54 


1.28 


0.08 


5 




15.086 


10 


43 


12.8952 


+11 


52 


35.685 


48 


78 


0.31 


0.73 


0.08 


6 


J104320. 5+114301 


9.1344 


10 


43 


20.5287 


+11 


43 


01.933 


10 


16 


1.48 


5.00 


0.08 


7 


J104328.4+113723 


8.6480 


10 


43 


28.4280 


+11 


37 


23.613 


5 


10 


2.43 


5.53 


0.08 


8 




5.2228 


10 


43 


36.7801 


+11 


43 


15.556 


21 


34 


0.72 


1.65 


0.08 


9 


J104342. 7+115209 


10.575 


10 


43 


42.7584 


+11 


52 


09.373 


2 


4 


4.33 


11.82 


0.08 


10 




2.8629 


10 


43 


47.3236 


+11 


40 


54.775 


35 


57 


0.42 


0.81 


0.08 


11 


J104349. 1+113628 


6.1266 


10 


43 


49.1642 


+11 


36 


28.452 


3 


5 


2.74 


4.00 


0.09 


12 




1.9497 


10 


43 


49.8033 


+11 


41 


58.791 


43 


70 


0.43 


0.85 


0.08 


13 




2.8251 


10 


43 


50.6611 


+11 


39 


59.677 


35 


56 


0.42 


0.78 


0.08 


14(R1) 


J104355. 2+114129 


0.9541 


10 


43 


55.2618 


+11 


41 


29.337 


2 


4 


4.11 


4.08 


0.07 


15(R2) 


J104359. 5+114332 


1.3864 


10 


43 


59.5464 


+11 


43 


32.638 


3 


5 


3.53 


22.04 


0.07 


16(R3) 


J104401.4+114323 


1.4793 


10 


44 


01.4482 


+11 


43 


23.623 


21 


18 


1.38 


15.93 


0.06 


17 




2.9359 


10 


44 


08.4135 


+11 


40 


54.837 


38 


58 


0.39 


1.09 


0.08 


18 




3.0888 


10 


44 


10.3120 


+11 


42 


08.643 


35 


57 


0.42 


0.50 


0.08 


19 


J104417.4+114614 


6.2824 


10 


44 


17.4533 


+11 


46 


14.678 


1 


1 


62.77 


125.40 


0.21 


20 




6.6425 


10 


44 


17.5254 


+11 


46 


46.157 


15 


24 


0.96 


1.79 


0.09 


21 


J104423.4+114458 


6.8677 


10 


44 


23.4286 


+11 


44 


58.340 


3 


3 


14.05 


32.41 


0.32 


22 


J104429. 2+113811 


8.7224 


10 


44 


29.2655 


+11 


38 


11.349 


7 


6 


4.05 


13.05 


0.14 


23 


J104430.6+113811 


9.0252 


10 


44 


30.6493 


+11 


38 


11.184 


1 


1 


21.12 


88.18 


0.13 



used a 4s integration time. After flagging bad data, the to- 
tal, usable time on NGC3351 was 15.4 hours, which is ex- 
pected to give a noise level of cr rm s = 53 /iJy beam -1 under 
good conditions. The point-like QSO OQ208 was used as a 
bandpass calibration source and to set the flux scale on all 
baselines. Its flux density was found to be 0.955 =b 0.008 Jy 
by comparison with 3C286 on the shorter baselines, where 
the latter is unresolved and has a flux density of 14.73 Jy at 
1.42 GHz (Baars et al. 1977). 

We followed the standard MERLIN continuum data re- 
duction procedure (MERLIN User Guide, Diamond et al. 
2003 1 ), using the AIPS package 2 for calibration and imag- 
ing. We applied the flux scale, bandpass and time- variable 
phase and amplitude solutions from the reference sources 
to NGC3351 and split out the data. We retained 15 1- 
MHz channels. This spectral resolution gives the most strin- 
gent restriction on the field of view, of about 10 percent 
smearing at 140 arcsec radius. We used separate facets of 
a few arcmin in all mappings in order to avoid distortions 
due to sky curvature and natural weighting, giving a beam 
size of 0.29 x 0.17 arcsec with a position angle of 23°. We 
searched for confusing sources by making a grid of 1681 
100-arcsec maps, without cleaning, over a 1-degree square 
field. We then re-imaged and cleaned all fields containing 
sources brighter than 2 mJy, using the clean components as 
a model for self-calibration of the whole dataset. We car- 
ried out further rounds of imaging and self-calibration, test- 
ing different strategies such as establishing that restricting 

1 http:/ /www. merlin. ac.uk/user_guide/OnlineMUG/ 

2 http://www.aips.nrao.edu/ 



self-calibration to the inner sources gave poorer results. We 
established that sidelobes fell to < 5 percent of the peak 
at distances > 0.5 arcmin from a source and identified all 
sources brighter than 1 mJy within the central 12 arcmin. 
Fainter sources should not produce sidelobes above the ex- 
pected arms in the central region. 



3 RESULTS 

Table 1 gives the properties of the 23 radio-continuum 
sources at 1.42 GHz found by MERLIN within a 35 x 22 ar- 
cmin region around NGC3351. The first column identifies 
the sources. The three objects that are closest to the centre 
(14, 15, and 16) will be referred to as Rl, R2, and R3 in what 
follows. Fifteen of the MERLIN sources have a counterpart 
in the VLA FIRST (Very Large Array, Faint Images of the 
Radio Sky at Twenty-cm) survey catalogue (White et al. 
1997) within 1 arcsec radius; names of the FIRST counter- 
parts are quoted in column 2. Column 3 shows the radial 
offset from the MERLIN pointing centre, and columns 4 to 
7 give the coordinates of each source and their correspond- 
ing errors. Peak surface brightness, integrated flux densities, 
and local noise levels cr rms are listed in columns 8, 9, and 
10, respectively. The distribution of the MERLIN sources is 
shown in Fig. 1 overplotted on an Ha image acquired with 
the Palomar 48-inch Schmidt telescope and obtained from 
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Figure 1. Distribution of the MERLIN sources in a 35 x 22 arcmin 2 field centred at 0:2000 = 10 h 43 m 57!7, 5 2 ooo = +ll°42 / 14 / /0 su- 
perimposed in an Ha image from Palomar 48-inch Schmidt telescope. Crosses and small circles indicate the MERLIN and FIRST 
sources, respectively. The MERLIN sources are labelled according to Table 1. The big circle indicates the central 100 arcsec (in radius) 
of NGC3351, and the ellipse is the D25 isophote of this galaxy as defined by de Vaucouleurs et al. (1991). The orientation is north up, 
east to the left. [See the electronic edition of the Journal for a colour version of this figure.] 



the Digitized Sky Survey 3 (DSS) through the NASA/IPAC 
Extragalactic Database 4 (NED). 

We found six sources inside the D25 isophote whose in- 
tensity peaks are brighter than the local 6 cr rm s confidence 
limit of 310 — 340 /iJybeam -1 , corresponding to a bright- 
ness temperature of ~1000K. Outside this ellipse (D25), 
the noise increases, reaching about 0.07-0.08 mJy beam -1 
in the regions out to ^15-20 arcmin radius which are af- 
fected by bandwidth and other smearing, making accurate 
deconvolution difficult. We visually examined each candi- 
date source, rejecting those which seemed spurious. Arte- 
facts did not exceed 6 cr rms within the primary beam FWHM, 
except for clearly-recognisable sidelobes reaching 10a rms 
within ~ 1 arcmin of some of the brightest outlying sources 
In the outer part of the field, we retained fifteen sources ex- 
ceeding 10 arms (~ 700 — 800 mJy beam -1 ) and two fainter 
sources (5 and 18) which are conspicuously isolated from any 
contaminating sidelobes. The distribution of radio sources 
close to the centre of NGC3351 is shown with crosses in 
Figure 2. This false colour image is made using two Spitzer 
images at 8.0 and 3.6 /xm (red and green channels, respec- 
tively) and a continuum-subtracted Ha image (blue channel) 



3 http:/ /stdatu. stsci.edu/cgi-bin/dss_form 

4 http: / / nedwww.ipac.caltech.edu/ 



from CTIO 4.0 m Mosaic2 (correction to the celestial coordi- 
nates of this image was provided by Douglas Swartz, private 
communication) . 

MERLIN detected 15 (about 63%) of the 24 sources 
listed in the FIRST survey catalogue within the field of view 
shown in Fig. 1. The area of the MERLIN beam is < 1% 
of the FIRST resolution of ~ 5 arcsec, giving a brightness 
sensitivity for MERLIN > 1000 K compared with 20-30 K 
for FIRST. In addition, the shortest MERLIN baselines cor- 
respond to a maximum imageable angular scale of ~ 2 arc- 
sec. All the non-detected sources are fainter than 7 mJy per 
FIRST beam. On the other hand, we have discovered 8 new 
sources that were not previously catalogued in FIRST, due 
to our greater sensitivity to compact emission. The posi- 
tion errors quoted in Table 1 are the noise-based relative er- 
rors of the source peaks. Astrometric uncertainties also arise 
from the assumed positions of the phase reference source 
(~ 3 mas) and the MERLIN antennas (~ 10 mas), as well as 
from the phase reference-target separation (~ 5 mas) , giving 
a total astrometric uncertainty of about 12 mas. The sources 
at increasing distances from the centre, over a few arcmin, 
may suffer from additional errors in the peak position due 
to smearing. 

The best images of the central region were achieved 
by using all sources identified in a 1-degree region for self- 
calibration and subtracting the clean components for sources 
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Figure 2. Superposition of the Spitzer images at 8.0 /im (red) 
and 3.6 /im (green) with H a (blue) emission distribution from 
NGC 3351. The yellow crosses indicate the radio sources projected 
over the galaxy, labelled according to Table 1. The big red circle 
indicates the central lOOarcsec (in radius) of NGC 3351, and the 
cyan ellipse is the D25 isophote of this galaxy as denned by de 
Vaucouleurs et al. (1991). The orientation is north up, east to the 
left. 



further than 2 arcmin from the centre. After all calibra- 
tion and subtraction, we made a cleaned image of the re- 
gion within 100 arcsec of the centre. Using natural weighting 
gave cr rms = 57 fiJy beam -1 . We tried giving the longer base- 
lines lower weight (tapering) and used a restoring beam of 
0.5 x 0.5 arcsec, giving cr rms = 52 /xJy. Although these levels 
are close to the predicted noise (allowing for the low declina- 
tion of NGC 3351), even the best images contained ripples of 
magnitude ~ 6 cr rm s due to the severe distortions introduced 
by bright sources at large offsets (which have peaks of up to 
a few Jy in FIRST) that could not be adequately subtracted 
or cleaned out. 

Three sources, abbreviated by us as Rl, R2 and R3, 
were found within the central 100 arcsec in radius (see Fig- 
ures 2 and 3) . All of them are brighter than 1 m Jy, and 
both R2 and R3 are significantly extended. Rl seems to be 
very compact, and MERLIN measured a higher flux density 
(4.1 mjy) than FIRST (1.9 mjy), suggesting not only that 
this source is smaller than the MERLIN beam, but also that 
it is compact enough to be highly variable, having approxi- 
mately doubled in flux between the date of the FIRST sur- 
vey and that of the MERLIN observations. The FITS header 
of the observations in the VLA archive shows that they were 
made in Dec 1999 and Jan 2000. That means that the size 
of Rl is constrained to be about 2.4 pc, corresponding to a 
light-travel diameter of 8 yr. 

As can be appreciated in Figure 3, these sources display 
several faint, yet observable, substructures. Rl displays a rel- 
atively compact main knot, and the other features are too 
faint to be sure they are real. R2 has an unresolved knot to 
the southwest, not necessarily compact, but its components 
are too faint and irregular, or too blended, to measure their 
extent reliably. R3 shows three fainter knots. Although the 
intensity peak (in mjy per beam) of the brightest one is al- 



most equal to that of the main knot, its integrated flux is 
less than a half. Table 2 lists the coordinates, peak surface 
brightness, integrated flux densities, and sizes for each knot, 
compared to the analogous values given in the FIRST survey 
(White et al. 1997). There are a few other features extended 
along the beam angle, but they seem to be artefacts and thus 
we have not listed them in the table. All fluxes and sizes, 
as well as their corresponding errors, have been estimated 
from the 4, 6 and 8cr rms level contours. Sizes are given af- 
ter deconvolution of the restoring beam, and their errors are 
smaller than 0.02 arcsec (5 in position angle) for sources 
brighter than 1 mjy. For fainter sources, the errors increase 
to 0.05 arcsec and 20 for sizes and position angles, respec- 
tively. The discrepancy between the two quoted position an- 
gles of R2 between the FIRST and MERLIN measurements 
(see middle panel of Figure 3) is partly due to the fact that 
the MERLIN PA is for the brighter peak only, whilst the 
FIRST PA is derived from a fit to the entire patch of emis- 
sion. The position angle between the two MERLIN peaks is 
136 deg, which agrees with the FIRST PA within the errors. 
FIRST does not give individual position angle errors. We 
have downloaded the image from the FIRST archive 5 and 
fitting R2 and R3 we have obtained an uncertainty of 14 and 
18 deg, respectively. 

Rl is the only source that has a clear counterpart at 
other wavelengths, from X-rays to radio. It has already 
been observed and identified as an X-ray point source us- 
ing Chandra data (Sll in the catalogue of Swartz et al. 
2006). Galex at near- and far-UV, as well as Spitzer-MIPS 
at 70 and 160 /im in the far-IR, do not show a clear coun- 
terpart due to their lower spatial resolution, leading to con- 
fusion and blending with nearby objects. Images of the re- 
gion around Rl in different bands are shown in Figure 4: 
X-rays from the Chandra Data Archive 6 , Ha and IR images 
from NED, HST-WFPC2 images from the Hubble Legacy 
Archive 7 (HLA) and the Multimission Archive at STScI 8 
(MAST). The X-ray image has been built by combining the 
three ~ 40 ks images available in the archive using the CIAO 
3.4 software (Fruscione et al. 2006). For R2 and R3, we find 
no obvious counterpart at any wavelength (see Figure 5). 
In fact, they seem to be located in regions devoid of emis- 
sion in other bands. There is no correlation between these 
MERLIN radio sources and other catalogues, such as the at- 
las of Hn regions by Hodge & Kennicutt (1983) or Bradley 
et al. (2006). A general search in NED for entries around the 
positions of R2 and R3 also yielded negative results. They 
do not have satisfactory cross-identifications in SPECFIND 
(Vollmer et al. 2010), although either or both could be as- 
sociated with 611 and 4850 MHz sources observed at lower 
resolution. The high flux densities observed at 611 MHz hint 
at these sources having a non-thermal origin. 

We found no further sources brighter than ~ 6 cr rms in 
the central region. In particular, the two nuclear sources 
catalogued in FIRST are resolved out in the present data. At 
larger radii, MERLIN observed 3 additional sources above 
the 6 cr rms limit that are projected inside the D25 isophote of 



5 http://sundog.stsci.edu/ 

6 http: / /cda. harvard.edu / chaser / mainEntry.do 

7 http:/ /hla. stsci.edu/hlaview. html 

8 http://archive.stsci.edu/ 
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Figure 3. Radio continuum MERLIN map of the sources Rl, R2 and R3 (left, middle and right panel, respectively). The scale is the 
same in all the panels. Level contours correspond to -1, 1, 2, 4... x 0.2 mjy beam -1 . The crosses represent the major and minor axes of 
the deconvolved components fitted to the MERLIN images and the ellipses correspond to the FIRST deconvolved source shapes (the 
dashed circle represents the unresolved source). The orientation is north up, east to the left. [See the electronic edition of the Journal 
for a colour version of this figure.] 



Table 2. Radio continuum parameters at 1.42 GHz from the extended sources. The first column gives the adopted reference name. From 
column 2 to 6 we list the peak position, the peak intensity, integrated flux density and sizes (major axis x minor axis, position angle) 
of the different knots of these sources. The last two columns give the values listed in the FIRST survey (White et al. 1997). 



Region 






MERLIN 








FIRST 




R.A. 


Dec. 


Peak 


Int. 


Size 


Int 


Size 




(J2000) 


(mJy beam -1 ) 


(mJy) 


(arcsec x arcsec,deg. ) 


(mJy) 


(arcsec x arcsec, deg.) 


Rl 


10 43 55.2618 


+11 41 29.337 


4.11 


4.08 + 0.06 


(compact) 


1.91 


(compact) 


R2 


10 43 59.5480 


+11 43 32.572 


3.82 


23.76 + 0.37 


0.65 x 0.35,101 


22.12 


1.64 x 0.93,151.3 




10 43 59.5338 


+11 43 32.810 


3.53 


3.53 + 0.06 


(unresolved) 






R3 


10 44 01.4465 


+11 43 23.646 


1.50 


10.57 + 0.41 


0.85 x 0.28,103 


23.19 


2.19 x 1.81,109.6 




10 44 01.5560 


+11 43 23.228 


1.46 


4.93 + 0.23 


0.39 x 0.30,17 








10 44 01.4783 


+11 43 23.482 


0.92 


3.84 + 0.27 


0.51 x 0.26,144 








10 44 01.3798 


+11 43 23.832 


0.47 


1.29 + 0.19 


0.32 x 0.24,152 







NGC3351. These sources (entries 10, 12, and 13 in Table 1, 
located at 2-3 arcmin from the centre of the observed field) 
were not detected in the FIRST survey nor in the VLA 
observations reported in Paladino et al. (2006). 



4 DISCUSSION 

Early- type spiral galaxies such as NGC 3351 provide an ideal 
laboratory to study the very different environments in which 
star formation may take place, from the peculiar conditions 
of the CNSFRs (e.g. Hagele et al. 2009, 2010, and refer- 
ences therein) to the outskirts of the galaxy disk. Maps of 
radio-continuum emission can reveal the presence of SNRs 
and Hll regions through their synchrotron and thermal emis- 
sion, respectively. These are closely related to the evolution 
of massive stars and are, therefore, excellent tracers of the 
recent star formation activity within the galaxy (see e.g. 
Muxlow et al. 1994). Moreover, young massive star- forming 
regions have very characteristic signatures in the radio band, 
providing strong constraints on the age of the most recent 
starbursts (see e.g. Rosa-Gonzalez et al. 2007). 

We have divided the galaxy in several regions, corre- 
sponding to different star formation regimes: the centre, 
which we define as the central 600 pc (12 arcsec radius) 
features the highest densities and metallicities, Z ~ Z . 



The CNSFR complexes identified in Planesas et al. (1997) 
have been considered separately from the diffuse nuclear star 
formation. Finally, the galaxy disk (r > 600 pc, r > 12 arc- 
sec) represents a more typical, low-density environment with 
low to moderate metallicities. The main properties (diame- 
ter, Ha luminosity, electron density and temperature, filling 
factor, thermal optical depth and flux density, star forma- 
tion rate, total number of SNRs, expected number above our 
detection threshold, and total synchrotron flux density) of 
each region are summarised in Table 3. The eight CNSFRs 
are listed in the same order as in Planesas et al. (1997). 

Diameters (D) of the CNSFR complexes are estimated 
from the HST F606W image, and Ha luminosities [L(Ha)] 
are taken from Planesas et al. (1997), corrected for inter- 
nal extinction using the colour excess [E(B-V)] estimated by 
Perez-Olea (1996) and assuming the Galactic extinction law 
of Miller & Mathews (1972) with R v = 3.2 (Hagele 2008). For 
the rest of the nuclear region (excluding the CNSFRs), the 
star formation rate derived from the Br7 luminosity (Pux- 
ley et al. 1990) has been converted to an Ha luminosity 
by means of the empirical relation given by Rosa-Gonzalez 
et al. (2002). For the galactic disk, the SFR measured by 
Leroy et al. (2008) has been used. 

Electron densities (n e ) of the emitting gas in the CNS- 
FRs have been computed from the [Sn]AA6717/6731 A line 
ratio, and electron temperatures refer to T e ([Oin]) (Diaz 
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Figure 4. Images of 1.15' x 0.75 ; around Rl. From top, left to right: Chandra 0.1-8 keV, HST-WFPC2 F439W (wide B), HST-WFPC2 
F606W (wide V), CTIO Ha, HST-WFPC2 F814W (wide I), Spitzer-IRAC 4.5 /im, Spitzer-IRAC 8 /xm, Spitzer-MIPS 24 /xm. In all panels 
crosses show the two detected knots of Rl, diamonds represent the x-rays sources detected by Chandra with the nomenclature defined 
by Swartz et al. (2006). In all the cases the orientation is north up, east to the left. [See the electronic edition of the Journal for a colour 
version of this figure.] 
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Figure 5. Images of 1.15' x 0.75' around R2 and R3. From top, left to right: Chandra 0.1-8 keV, HST-WFPC2 F450W (wide B), HST- 
WFPC2 F606W (wide V), CTIO Ha, Spitzer-IRAC 3.6 /mi, Spitzer-IRAC 5.8 /mi, Spitzer-IRAC 8 /mi, Spitzer-MIPS 24 /mi. In all panels 
crosses show the locations of R2 and R3, diamonds represent the x-rays sources detected by Chandra with the nomenclature defined by 
Swartz et al. (2006). In all the cases the orientation is north up, east to the left. [See the electronic edition of the Journal for a colour 
version of this figure.] 
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Table 3. Main properties of each region: diameter, Ha luminosity, electron density and temperature, filling factor, thermal optical depth 
and flux density, star formation rate, total number of SNRs, expected number above our detection threshold, and total synchrotron flux 
density. 



Region 


D 

(arcsec) 


L(Ha) 
(10 39 ergs" 1 ) 


n e 
(cm- 3 ) 


T e 
(10 4 K) 


/ 


-brems 
T 1.42 GHz 


cbrems 
°1A2 GHz 

(mJy) 


SFR 
(M yr- 1 ) 


NsNR 


,.(0.3 mJy) 
SNR 


qrsynchr 
^1.42 GHz 

(mJy) 


CNSFRs 


1.00 


5.43 


360 


0.48 


0.021 


0.212 


0.30 


0.060 


4.4 


0.119 


8.38 




1.00 


2.83 


440 


0.45 


0.007 


0.165 


0.16 


0.031 


2.3 


0.062 


4.37 




1.00 


7.00 


430 


0.45 


0.018 


0.296 


0.36 


0.077 


5.7 


0.154 


10.80 




0.80 


2.32 


310 


0.44 


0.022 


0.145 


0.13 


0.026 


1.9 


0.051 


3.58 




0.60 


1.03 


360 


0.54 


0.021 


0.106 


0.06 


0.011 


0.8 


0.023 


1.59 




1.50 


0.75 


360 


0.55 


0.001 


0.034 


0.05 


0.008 


0.6 


0.017 


1.16 




1.00 


1.80 


410 


0.46 


0.005 


0.114 


0.10 


0.020 


1.5 


0.040 


2.78 




0.60 


0.54 


400 


0.50 


0.008 


0.079 


0.03 


0.006 


0.4 


0.012 


0.83 


Centre 


24.00 


23.60 


10 


1.00 


0.016 


0.001 


2.14 


0.260 


19.2 


0.519 


36.41 


Disk 


400.00 


45.50 


1 


1.00 


0.001 


0.000 


4.13 


0.500 


125.1 


1.001 


70.20 



et al. 2007). For the eighth region, values of 400 cm -3 and 
5000 K have been assumed. The ambient ISM densities of 
the central diffuse region and the disk are 10 and 1 cm -3 , 
respectively, with a temperature of 10 4 K. The associated 
filling factors (/) are given by expression (A5). As explained 
in Appendix A, the predicted optical depth (t^qhJ and 
flux density (S'i^ghz) at 1.42 GHz arising from thermal 
bremsstrahlung take into account self- absorption, although 
the effect is found to be mild at this frequency. 

Star formation rates have been estimated from the Ha 
luminosity (Rosa-Gonzalez et al. 2002), and the expected 
number of SNRs (A/snr) has been computed according to 
equation (B4), assuming an ambient ISM density of 10 cm -3 
for the CNSFRs. The number of SNRs brighter than 0.3 mJy 
(^snr mJy ' > ) is given by expression (B10), and the total syn- 
chrotron emission from the whole population of unresolved 
supernovae (^J^ghz) nas been estimated from the empiri- 
cal relation between the star formation rate and the 1.4 GHz 
luminosity density reported by Yun et al. (2001). 

As shown in the previous section, our analysis of the 
sub-arcsecond spatial resolution radio map of NGC3351 
shows the presence of 23 sources in the L-band. Three of 
them (Rl, R2 and R3) are located within the central 100 arc- 
sec in radius. They have clear counterparts in the FIRST 
survey (White et al. 1997) and the VLA data analysed by 
Paladino et al. (2006), but only Rl seems to be visible at 
other wavelengths. It is unresolved by MERLIN, putting an 
upper limit of ~ 10 pc on its size, but its variability (Sec- 
tion 3) suggests that it is < 3 pc. As shown in Figure 5, 
neither R2 nor R3 display significant emission in X-ray, UV, 
optical, or IR images. These sources are significantly ex- 
tended (see Figure 3). They are separated by about half an 
arcmin and could possibly be associated with jets from a 
background AGN. FIRST contains about 7 sources deg -2 
brighter than 20 mJy; the probability of finding one of them 
within a region of 100 arcsec radius chosen at random is 
only ~ 2 per cent, but their observed fluxes and extents 
are difficult to reconcile with a stellar origin. Moreover, 
the chance of two unrelated sources of such high bright- 
ness occurring so close together is extremely low, favouring 
the background AGN scenario. Thus, the only radio point 
source above 0.3 mJy that may be related to star formation 
in NGC3351 is Rl. 

In particular, we have not detected any point source 
in the nuclear region. The thermal emission predicted for 



each of the CNSFR complexes ranges from 0.03 to 0.36 mJy, 
slightly below our detection limit 9 . The estimated fluxes are 
fairly robust with respect to the adopted sizes. Smaller diam- 
eters, appropriate for the individual knots (see e.g. Hagele 
et al. 2007), yield higher filling factors and optical depths, 
smaller areas, and, then, similar fluxes. On the other hand, 
a number of supernovae is expected to be present in these 
regions. The expected probability of finding a SNR brighter 
than 0.3 mJy in the circumnuclear ring is relatively high 
(~ 50 per cent), and the synchrotron flux associated to 
the unresolved population within each CNSFR (last column 
of Table 3) would have been easily detected by MERLIN. 
Therefore, the observed lack of emission suggests that these 
systems are younger than the lifetimes of the most massive 
stars (see e.g. Rosa-Gonzalez et al. 2007). This seems also 
to be the case in similar CNSFRs observed in other galaxies 
(e.g. Dors et al. 2008, in NGC6951). 

Regarding the star formation activity in the rest of the 
central region, the predicted number of bright SNR is also 
~ 0.5, and the total synchrotron flux is estimated to be 
about 36 mJy. Indeed, VLA observations at a minimum 
resolution of 4.5 x 4.1 arcsec reveal a large central core, 
about 20 arcsec in size, featuring a double emission peak 
with maximum surface brightness about 10 10 mJy sr _1 (Pal- 
adino et al. 2006), ~ 0.05 mJy per MERLIN beam (close to 
the rms noise of our observations) if the emission was uni- 
formly distributed. If star formation takes place in compact 
knots, similar in size to those in the CNSFRs, this emis- 
sion ought to display conspicuous small-scale structure that 
should be detected by MERLIN. However, MERLIN is only 
sensitive to scales smaller than ~ 2 arcsec, which is ~ 100 
pc at the distance of NGC3351. Therefore, the relativistic 
electrons generated by the supernova explosions must have 
diffused over distances longer than our maximum resolution 
radius of ~ 50 pc, and the true diffusion scale would thus be 
somewhere in between this value and the 300 pc (6 arcsec) 
where the radio-IR correlation is observed by Paladino et al. 
(2006). 

The total flux density at 1.42 GHz, integrated over the 
whole galaxy, is predicted to be ~ 0.14 Jy. According to 
equation (B10), the expected number of radio sources above 



9 For the requested observing time of 48 hours, the expected 
noise level would have been cr r ms = 0.03 mJy under good weather 
conditions. 



10 G. F. Hdgele et al. 



0.3 mJy in the disk of NGC 3351 is of the order of one, con- 
sistent with the interpretation that Rl is a young SNR. 
However, an increase in flux of a factor of two between the 
FIRST observations and our data seems difficult to justify 
in this scenario. Two possible explanations would be that 
FIRST caught the supernova near the time of the explosion, 
where the radio continuum was still rising, and we are now 
observing a very young remnant, as suggested by its high lu- 
minosity. An alternative possibility would be that the radio 
flux has increased due to density variations in the ambient 
medium. Some Galactic nebulae, such as Cas A (Erickson 
& Perley 1975; Read 1977; Chevalier et al. 1978) and GK 
Per (Anupama & Kantharia 2005) also show anomalous in- 
creases in flux density, albeit at lower frequencies. Ideally 
SNRs could be identified based on the radio spectral index, 
and future observations will be sought based on the results 
of the current work. 
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5 SUMMARY 

We present high spatial resolution radio continuum obser- 
vations at 1.42 GHz acquired with the Multi-Element Radio 
Linked Interferometer Network (MERLIN) in and around 
NGC 3351. We found 23 radio sources in a field of 35 by 
22 arcmin centred in this galaxy. 6 of them are projected 
within its D25 isophote, and 3 are located inside the cen- 
tral lOOarcsec in radius. Two of these three are significantly 
extended, and the other one is the only source with a previ- 
ously detected counterpart at other wavelengths. Our results 
suggest that this radio source might correspond to a young 
supernova remnant, while the other two are probably related 
to jets from a background AGN. 

We do not detect individual supernovae or SNRs in the 
central region of the galaxy. The CNSFRs are too young 
(less than a few Myr) to host supernovae, and the diffusion 
length of the relativistic electrons in the ISM seems to be 
larger than our maximum resolution of 50 pc. Detecting the 
thermal bremsstrahlung emission from the circumnuclear 
Hn regions requires deeper observations. e-MERLIN will al- 
low NGC 3351 to be mapped at a sensitivity of cr rms ~ 10/iJy 
beam -1 in a similar observation, thanks not only to the in- 
creased sensitivity provided by 0.5 MHz bandwidth around 
1.5 GHz but also to the new correlator which will allow the 
entire primary beam to be imaged without distortion. 
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APPENDIX A: THERMAL 
BREMSSTRAHLUNG FROM Hn REGIONS 

Thermal bremsstrahlung from Hn regions can be a potential 
signature of star formation in the radio band. The optical 
depth through an ionised hydrogen cloud can be approxi- 
mated by (Mezger & Henderson 1967) 

where the emission measure is defined as the integral along 
the line of sight 

EM = J n e n { dltt Jul til (A2) 

with n e and rti denoting electron and ion density, respec- 
tively. For a pure hydrogen, fully ionised plasma, n e = rii. 
Including helium would introduce a correction of the order 
of 10 per cent, and the contribution of heavier elements is 
negligible. 

The overall extent of the CNSFRs in NGC 3351 is 
~ 100 pc, but most of the emission comes from a few com- 
pact knots of (1 — 5) pc radius. The characteristic electron 
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density of the emitting gas is accurately measured by the 
ratio between the optical sulphur lines [Sn] A A671 7/6731 A, 
with typical values of ~ 400 cm -3 (Diaz et al. 2007). The 
volume filling factor / of the (unresolved) knots can be de- 
rived from Ha luminosity of the entire region 



L(Hq) = E a a a n 2 e f^D 3 



(A3) 



where E a = 3 x 10 12 erg is the energy of the Ha transition, 



1.17 x 10" 



(A4) 



V10 4 K J 

is the rate coefficient, and a cylindrical geometry with both 
diameter and depth D has been assumed for simplicity. The 
uncertainty associated with this assumption is high (e.g. the 
difference with respect to a perfect sphere of diameter D is 
a factor of 1.5), and it constitutes an important source of 
error in the following estimates. With this caveat in mind, 
the filling factor of the emitting gas is approximately given 
by the expression 

L(Ha) _^ (Ji^X 2 (D^ :i ' " ^"• S7 
^cm -3 / \pc 



/ ~1.2 



(A5) 



^lO^ergs- 1 J Vcm- 3 ; \pc J V 1Q4K , 

If there is a small number N of identical knots within the 
region, so that they do not overlap along the line of sight, 
the total emission measure would be 

,1/3 



EM - 



D 



and the total area that they subtend on the sky is 

/X 2/3 



N- 



N 



(A6) 



(A7) 



whereas for large N, assuming that the knots are randomly 
distributed, 



EM = nlfD (A8) 
and 

A = \D 2 (A9) 

The two cases N = 1 and N — >• oo can be considered as 
reasonable upper and lower limits to the emission measure, 
respectively. The radio optical depth is thus bracketed by 

01 / ^ \-1.35 

r„ « 3.3xl0- 7 (— ) 

x2 / D 



10 4 K 



X 



(— V 

Vcm -3 / 



pcj 



(A10) 



with 1/3 < rj < 1. Taking the geometric average, 77 = 2/3, 
the relative error incurred in the optical depth is smaller 
than a factor f 1 ^ 3 . 

The radio-continuum surface brightness is given by the 
expression 

E„ = B v (1 - e _T ") (All) 
where B u w 2kT ^ [ s the Planck function for a black body 



at temperature T in the Ray leigh- Jeans approximation, and 
the flux density at the Earth is 

Jy 



d 2 



2.4 x m 



VGHz/ \10*K 



D 



I 1 ' 11 (1 - e- T -()A12) 



In the optically thin case, Y, u « B u r u implies 



— « 79 l 
Jy VGHz, 



Tie \ 2 fD 



Vcm A J 



V10 4 K/ \ct y vein 
and one recovers the well-known relation 



pc 



^«l.lxl0 9 

Jy 



(ghz) 



10 4 K 



.F(Ha) 



(A13) 



(A14) 



between the radio continuum and Ha fluxes. In the general 
case, the radio-continuum flux density may be estimated by 
combining expressions (A5), (A10), and (A12). 



APPENDIX B: SYNCHROTRON EMISSION 
FROM SUPERNOVA REMNANTS 

Assuming that all stars more massive than Msn -8M0 ex- 
plode as core collapse (Type II) supernovae, the number 
of such events per unit time is determined by the SFR M* 
and the initial mass function (IMF) of the stellar population, 



^sn(£) 



J^mfPim) dm 



(Bl) 



where t* (m) is the lifetime of a star of mass m. For a Salpeter 
(1955) IMF, ip(m) oc m -2-35 , with minimum and maximum 
stellar masses M m in 00 0.1 M© and M max ~ 100 M©, respec- 
tively. If the SFR is constant over the period t*(Msn) ~ 
40 Myr that the least massive supernovae take to explode, 



^sn(£) 



f A ^ max m ijj(m) dm 

iw min 



7.4 xlO- 3 ^^ (B2) 



SNRs can be detected as compact synchrotron sources 
during the characteristic time 

4/17 



-9/17 



^-3.4xl0 4 (-f^) ' (-^V ' (B3) 
yr V10 51 erg / Vcm- 3 / v J 

of the adiabatic expansion phase, where Esn denotes the 
energy input by the original supernova explosion, and n is 
the density of the ambient gas (Woltjer 1972). The total 
number of SNR observable at a given time in the radio band 
would be 

A^SNR(t) = / ^sn(£ ; ) dt' « ^sn(£)t S nr 

J t — TSIVR 

/ 7? x-9/17 

* (^) (B4) 



250 



t — T SNR 

M*(t) 



M© yr -1 Vcm - 

assuming Esn — 10 51 erg. The approximation of constant 
SFR for the last 40 Myr and a density n ~ 1 cm -3 are 
probably valid for most regions of the galaxy disk, where a 
moderate star formation surface density of the order of £* ~ 
(1 — 3) x 10 -3 M0 yr -1 kpc -2 is observed without prominent 
signatures of intense starbursts (see e.g. Leroy et al. 2008). 
The expected supernova surface density throughout the disk 
is thus 



-'SNR ; 



250 



0.5 kpc 



(B5) 



M yr- 1 

and the total number of SNR in a corona between 0.6 
and 10 kpc (12 and 200 arcsec) from the galactic centre 
would be about 125, corresponding to an integrated SFR 
of 0.5 M yr" 1 . 
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The inner 600 pc of NGC3351 are also forming stars 
at a rate of ~ 0.5 M yr _1 . A large fraction (about 
0.24 M© yr _1 , see Hagele et al. 2007) can be accounted for 
by the young starbursts of the CNSFRs, yielding a diffuse 
contribution of ~ 0.26 M yr _1 . Although equation (B4) 
predicts a similar number of supernovae in the central re- 
gion than in the whole disk, the density of the ISM is ar- 
guably higher. Using n — 10 cm -3 , the expected number 
of SNR becomes ~ 40. Concerning the individual CNSFRs, 
the observed values of M* ~ (0.01 — 0.06) Mq yr _1 suggest 
that most bursts should contain SNRs. However, one has to 
bear in mind that the star formation activity may have var- 
ied considerably over the last 40 Myr. In the most extreme 
case, there is the possibility that all bursts are younger than 
3 Myr, and therefore no supernova could have exploded yet. 
On the other hand, it is likely that the gas inside and around 
the CNSFRs is considerably denser than 10 cm -3 , leading 
to a shorter radio lifetime tsnr and a smaller number of 
radio-emitting SNR. 

Regardless of the number of sources actually present in 
the galaxy, they must be bright enough in order to be in- 
cluded in our sample. A relation between the radio surface 
brightness T, u and the diameter D of SNRs has long been 
established both theoretically (Shklovskii 1960) and obser- 
vationally (Clark & Caswell 1976). At 1 GHz, Urosevic et al. 
(2005) obtain 

SlGHz S.84 x 10- 16 f-V 3 ' 2 (B6) 



W Hz -1 m~ 2 sr _1 \pc J 

combining data from different nearby galaxies. To estimate 
the flux density of an extended source, one has to multiply 
its surface brightness by the solid angle subtended by the 
MERLIN beam, 

Si. 42 GHz _ ( 0FWHmYv rr ( D \ 3 2 , . 

7T Si. 4 GHz = 55 (B7) 



mjy beam -1 \ 2 J ' \pc 

where we have used #fwhm = 0.2 arcsec and a spectral index 
a — 0.5 to convert 1 GHz fluxes to 1.42 GHz. At the adopted 
distance of NGC 3351, an angle of 0.2 arcsec corresponds to a 
physical size of ~ 10 pc and a flux density of 35 fiJy beam -1 , 
below the noise level of the present data. Therefore, any SNR 
above our detection threshold must appear as an unresolved 
point source with flux 



/xJy \ a J \pc 

For a source as bright as Rl (4 mjy), the expected diam- 
eter is only 0.20 pc, for which the SNR would still be in 
the free expansion phase. Assuming an average expansion 
velocity of ~ 3000 km s -1 , the age of this system should be 
about 35 years. Substituting this number instead of tsnr in 
equation (B4), one obtains 

iVsNR(Si.42GHz^4 mjy) k 0.25 y^*^ (B9) 

Lowering the threshold to 300 /iJy, the diameter increases 
to 1.72 pc, and the expected number of sources becomes 

iV S NR(Si.42GHz^0.3 mjy) a 2 ( B1 °) 



